physica status solidi, 2 February 2008 



O 

o 
o 

(N 
> 

o 



Exciton dephasing in quantum dots: 
Coupling to LO phonons via excited 
states 

E. A. Muljarov 1 '*, R. Zimmermann 2 

1 General Physics Institute, Russian Academy of Sciences, Vavilova 38, Moscow 1 19991, Russia 

2 Institut fur Physik der Humboldt-Universitat zu Berlin, NewtonstraBe 15, 12489 Berlin, Germany 

Received XXXX, revised XXXX, accepted XXXX 
Published online XXXX 



PACS 78.67.Hc, 71.38.-k, 03.65.Yz, 71.35.-y 



T3 

o 
o 



' Con'esponding author: e-mail muljarov@gpi.ru 



We have found a novel mechanism of spectral broadening and dephasing in quantum dots (QDs) due to the coupling to 
longitudinal-optical (LO) phonons. In theory, this mechanism comes into play only if the complete manifold of exciton 
levels (including those in the wetting-layer continuum) is taken into account. We demonstrate this nontrivial dephasing 
in different types of QDs, using the exactly solvable quadratic coupling model, here generalized to an arbitrary number 
of excitonic states. 
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While acoustic phonons are known to produce a spec- 
tral broadening of the exciton zero phonon line (ZPL) in 
QD spectra, it is widely believed that the coupling to bulk 
dispersionless LO phonons forms everlasting mixed states 
(exciton polarons) showing no line broadening Q]|2|. This 
is indeed true if the model is restricted to a limited number 
of exciton states. However, including a sufficiently large 
number of states results in spectral broadening and dephas- 
ing. We have recently demonstrated [3 1 this nontrivial and 
quite unexpected mechanism of dephasing using the ex- 
actly solvable quadratic coupling model earlier developed 
for acoustic phonons |4|. This new mechanism is in clear 
contrast with the acoustic phonon-induced dephasing as 
the latter provides a finite ZPL width due to the phonon- 
mediated coupling already to the next exciton level. 

In this contribution, we derive the quadratic coupling 
model for a single level (exciton ground state) from the 
standard linear coupling between exciton levels in a QD 
and bulk LO phonons. Exciton excited states then con- 
tribute effectively through virtual transitions. We have gen- 
eralized the exact solution of the model with two levels (0 
to an arbitrary number of states J3| . Here we show ex- 
plicitly a diagonalization of the quadratic Hamiltonian and 
then calculate the dephasing in different types of QDs. 

We derive the model of quadratic coupling to LO phonons 
from the standard Hamiltonian with a linear exciton-phonon 



coupling: 
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where luq is the dispersionless LO phonon frequency, E n 
is the bare transition energy of a single-exciton state \n) in 

are matrix el- 



a QD, M nm (q) oc ? - 1 (n|e* t » r - - e"» r "| 
ements of the linear exciton-phonon Frohlich interaction, 
and % = 1. This problem has been solved exactly in case 
of a few exciton levels J2), but generally it is not solv- 
able for an arbitrary number of levels because of the level- 
nondiagonal coupling which is, at the same time, of key 
importance for the dephasing. To preserve this coupling, 
on the one hand, and to bring the problem to an exactly 
solvable form, on the other, we make the following unitary 
transformation 



H' = e s He- s = H+[S,H} + [S, [S, H}}/2 + . 
where the transformation operator S has the form: 
S=Y, \n)(m\J2Mnm(q) 
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This transformation eliminated nondiagonal terms in first 
order @), by mapping them into diagonal ones which are 
however quadratic in phonon operators. Higher powers of 
Mnm (q) are also generated in the transformed Hamiltoni- 
an. However, far from resonance where the quadratic cou- 
pling model is valid, they lead only to small quantitative 
corrections and thus can be safely neglected. 

For the ground exciton state (n = 1) the resulting quad- 
ratic Hamiltonian takes the form: 



Hi = E[ + cjo K a i + \ 

+ 5^Mu(q) (a q + a! q ) 
q 

-\ E ^(P' q ) ( a P + a - P )( fl q + a -q) 



(4) 
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is a renormalized (polaron shifted) energy of the exciton 
ground state, and 

E n - Ei 
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is the kernel of the quadratic coupling, here expressed in 
terms of new functions F n (q). Note that the Hamiltonian 
Eq.© contains also the state-diagonal linear coupling M\\. 

The absorption spectrum is calculated as a Fourier trans- 
form of the linear response on a delta-pulse excitation. 



Pit) = 6(t)(e lHot e- lHlt ), 



(6) 



where a finite-temperature expectation value is taken over 
the unperturbed phonon system 



Hq = uJq ^ ( a q a, 
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To do so we follow the method of Ref. (2) and diagonalize 
Hi directly in a few successive steps. 

First of all, we make a shift transformation: 



(7) 



(/q will be given later). This is only to eliminate the lin- 
ear coupling Mil, all other terms in Eq.© remaining the 
same, except for a new renormalization of the transition 
energy to E\ = E[ + E^ with an additional polaron shift 



(2) 

£p ol = — 5Z q Mii(q)/ q . The next step is a rotation of 
the phonon basis, in order to make the quadratic coupling 



diagonal. To do this we diagonalize the Hermitian matrix 
(E n — Ei > ojq is assumed for simplicity): 



q v 
and introduce new phonon operators: 

B„ = X;Ci,(q)&q 



with 



(8) 



(9) 



(10) 



making the transformation orthogonal, 



Our calculation of the optical polarization P(t) also needs 
to adapt the unshifted operators a q to the form of Eq.fO by 
making a similar transformation, 



A v = ^]C^(q)a q; 



(11) 



which preserves diagonality of the unperturbed Hamilto- 
nian due to negligible LO phonon dispersion: 



H = w Y, ( A l A » + o 



(12) 



For the perturbed Hamiltonian Hi, the diagonalization pro- 
cedure is finalized as follows 



Hi-Ei=J2 



w 



(13) 

where we have switched to another Bose operators 

A„ = a+B„ + a-Bl=a+A v + a-Al+ip v (14) 
with 

(15) 



a » = m — ' 

4>u=Y C M)U = \ir E a(q)M u (-q) . 

U!q — A11 v — 

q q 

Now the linear polarization obviously becomes a prod- 
uct of contributions from independent phonon subsystems 
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(labelled by the index v): 



P(t) = 9(t)e 
X 
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where e\ = (I + l/2)i2 v and e% = (k + l/2)u are, re- 
spectively, eigenenergies of the partial Hamiltonians H\ v 
and77 0l y: 

H 1 -E 1 =^2H lv , H lv \l\U>=e\\liU>, (17) 

V 

H = Y,H 0v , H 0v \k;v)=e v k \k;v), (18) 

V 

f3 = l/ksT, and 7 = 2cosh(/3u;o/2). The projections 
z l k {v) = (/c; v> of the perturbed phonon states (A^) 
into the bare ones {A„) are calculated recursively (below 
the index v is omitted for simplicity; k, I > 0): 
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where the second equation can be used to determine the 
start values for the recursion (I = 0), together with a nor- 
malization condition YlkLo \ z k\ 2 = 1- 

In case of two levels (only one excited state n = 2 
is taken into account), the matrix A nm becomes a scalar 
[see Eq.©] and thus only one phonon mode is renormal- 
ized. This phonon mode is bound to the QD and has a new 
frequency Q given by Eq.dTsb. Due to the LO phonon de- 
generacy, all the other phonon modes are decoupled and 
thus do not contribute to the spectrum. In fact, the remain- 
ing phonon modes can be always orthogonalized (with re- 
spect to the bound modes) without change of their fre- 
quency. According to Eq.([T6ll, the spectrum represents a 
set of discrete lines with a fine structure splitting given by 
the difference between u>o and Q, see Fig.l (a). Increasing 
the number of exciton levels participating in the coupling, 
the number of bound phonon modes (having renormalized 
frequencies, Eq.(fT5b. different from ujq) also increases and 
thus delta lines in the spectrum become more dense, see 
Fig.l (b) and (c). Finally, if we include all exciton levels, 
these lines merge and produce a continuous spectral broad- 
ening, Fig.l (d). 

This nontrivial mechanism of dephasing can be under- 
stood in terms of bound phonon states [6|. As shown in 
our recent paper 0, the bound-phonon frequencies £l v 
are condensed just below ujq. Such an accumulation of in- 
finitely many levels becomes a mathematical clue for the 
continuous spectral broadening. There are even some simi- 
larities to the model of quadratic coupling to acoustic pho- 
nons, in which taking into account only one excited state 
already results in infinitely many bound-phonon modes with 
renormalized frequencies and finally in a spectral broaden- 
ing (4|. However, the flat dispersion of LO phonons makes 




Figure 1 Calculated absorption spectra of the ground ex- 
citon state in a CdSe QD in the model of infinite har- 
monic potentials with the distance of 80meV (35meV) 
between electron (hole) levels coupled to LO phonons 
(cj = 24meV). 



the polarization decay essentially non-exponential, in con- 
trast with purely exponential long-time behavior caused by 
acoustic phonons JUG]. 

The absorption spectra in Fig.l are calculated for a 
QD with CdSe parameters, in the model of spherical har- 
monic potentials having only discrete levels. Now we an- 
alyze how the absorption of the QD changes if, instead of 
the unlimited parabolic confinement, a more realistic po- 
tential (parabola truncated at r = R) is used, so that the 
electron (and the hole) has both discrete and continuum 
states, see insets in Fig. 2. Below R = 5.5 nm [Fig.2(a,b)] 
only one discrete level exists, while for R = 5.5 nm and 
7nm [Fig. 2(c) and (d)] there are, respectively, one and 
three excited discrete states (shown by green horizontal 
lines). All other excited states (not shown) belong to the 
wetting layer (WL) continuum. States in the continuum 
are modelled in a large sphere with rigid walls. Numeri- 
cally, it is checked that for the given accuracy, there is no 
dependence on the number of continuum states included 
in the model and on the radius of the large sphere. In the 
worst case we have to take into account about 1500 densely 
spaced states in a lOOnm sphere. Below R = 7nm the 
spectra for truncated parabolas (black curved) are obvi- 
ously different from that of the unlimited parabola (red 
curves), see Fig. 2 (a) to (c). However, the linewidth remains 
practically the same. Finally, for R > 7 nm, the two spectra 
almost coincide, Fig. 2 (d). 
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Figure 2 Black curves: absorption calculated for different shapes of QD potentials shown in the insets along with the 
ground state wave function (shaded area) and discrete levels (horizontal lines) for the electron. Red curves: the spectrum 
calculated in the model of unlimited parabolic potentials, i.e. the same as in Fig.l (d). Zero of energy is the polaron shifted 
ground state transition energy; T = 150 K. 



In conclusion, a novel mechanism of spectral broaden- 
ing in QDs due to the standard linear exciton-LO phonon 
coupling is found and demonstrated by mapping it into ex- 
actly solvable quadratic coupling model. Generalizing this 
model to an arbitrary number of exciton excited states, the 
resulting Hamiltonian is diagonalized exactly and the op- 
tical polarization and absorption are calculated in terms 
of QD-bound phonon modes. It is shown that only taking 
into account the full excitonic spectrum, i.e. the complete 
manifold of excited states including both confined and the 
WL-continuum states, results in a LO phonon-induced ex- 
citon dephasing and spectral broadening. In simulation of 
the absorption in deep QDs with a few confined exciton 
states, the WL continuum can be safely replaced by a dis- 
crete spectrum, while in shallow QDs the continuum states 
have to be properly taken into account. 
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